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ABSTRACT 


An  underwater  explosion  forma  a gas  bubble  which  pulsates,  generates  a 
changing  pressure  field,  and  Induces  flexural  vibrations  In  a floating  or  submerged 
ship.  E<tuatlons  are  derived  which  couple  the  motion  In  each  normal  flexural  mode 
to  the  pulsation  cf  the  bubble.  Some  special  caBes  of  target  response  are  discussed. 
The  analysis  is  shown  to  be  applicable  If  the  duration  of  the  compressive  phase, 
the  shock  wave,  is  small  compared  with  the  pulsation  period  of  the  bubble  and 
with  the  flexural  period  of  the  target.  Equations  are  also  derived  which  give  the 
mode  shapes  and  frequencies  for  free  vibration  In  water  in  terms  of  these  modes 
and  frequencies  for  free  vibration  In  air. 


I 

INTRODUCTION 


When  an  explosive  charge  la  detonated  under  water,  the  detonation 
products  form  a gas  bubble  which  pulsates  and  migrates  In  the  water,  and  thus 
Induces  a changing  pressure  field.  These  pressures,  In  turn,  can  act  upon  a 
submerged  or  floating  target,  causing  bodily  motions  and  elastic  vibrations. 
There  already  exists  an  extensive  literature  which  describes  both  experimental 
and  theoretical  Investigations  Into  the  motions  of  the  gas  bubble,  but  without 
any  reference  to  Its  effect  upon  an  elastic  target.' 

Th.ls  report  presents  a theoretical  analysis  of  how  a slender  target, 
such  as  a ship  or  submarine,  reacts  to  the  pressures.  The  purpose  of  the  In- 
vestigation was  to  establish  the  physical  mechanism  responsible  for  the  flexu- 
ral motions  and  to  show  quantitatively  how  these  motions  depend  upon  the 
elastic  properties  of  the  target  and  upon  the  size,  position,  and  motion  of 
the  gas  bubble.  The  method  used  In  the  analysis  Is  to  separate  the  target 
reactions  and  water  motions  Into  normal  modes  of  motion  and  then  to  show  how 
each  mode  Is  coupled  to  the  bubble  pulsation.  The  effects  of  the  Initial 
shock  wave  are  at  first  omitted,  but  these  effects  are  Included  In  a later 
stage  of  the  analysis. 

Although  the  theory  may  seem  to  be  unduly  mathematl  a 1 divorced 
from  practical  problems,  the  applications  of  this  theory  are  significant  and 
will  be  discussed  In  subsequent  reports. 


N«*f>rrtn‘**»n  nr** 
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NORMAL  MODES  OP  THE  BODY 


We  shall  consider  only  small  vertical  motions  of  an  elongated  body 
for  which  the  Instantaneous  displacement  Is  uniform  for  all  points  In  the 
same  transverse  section.  Such  motions  are  generated  by  rigid  vertical  dis- 
placements, by  rigid  rotations  about  a horizontal  axis  perpendicular  to  the 
length,  and  also  by  combined  bending  and  shear  vibrations  In  planes  perpen- 
dicular to  the  length.  For  a body  which  Is  long  In  comparison  with  Its  width 
or  depth,  there  exists  a set  of  natural  modes  of  motion  of  this  type,  l.e., 
motions  which  It  may  have  when  moving  freely.  In  vacuum,  without  applied  ex- 
ternal forces.  These  motions  may  be  described  by  specifying  the  vertical 
displacement  ^(x,/)of  a longitudinal  axis  of  the  body  as  a function  of  the 
time  t and  the  distance  x from  one  end. 

^(x,  /)  « ^i(z)cos  dit  » =0,1,2-  * • [ 1 j 

The  functions  4>x  (x)  describe  the  mode  shapes.  The  Index  i , which  labels  the 
various  mode  functions,  may  conveniently  be  taken  as  equal  to  the  number  of 
nodes  of  ^(x)over  the  length  of  the  body.  The  frequencies  are  the  natural 
circular  frequencies  In  air.  These  mode  shapes  and  frequencies  completely 
specify  the  flexural  properties  of  the  body  and  will  be  taken  as  known  func- 
tions. In  practice,  they  may  be  determined  by  calculating  the  eigenfunctions 
and  eigenvalues  of  some  differential  equation  which  governs  the  flexure,  or 
perhaps  by  experiment. 

In  general  the  mode  functions  0,  are  orthogonal  over  the  length  of 
the  target,  with  m(x)the  mass  per  unit  length  as  a weighting  function,  so  that 

f \l>,  rndx  =»  A/td,,  t.j -0.1,2  •••  l2j 

Jo 


where  A/,  Is  defined  as  the  generalized  mass  associated  with  the  i‘*  mode. 
Also  these  mode  functions  form  a complete  set,  so  that  the  most  general  dy- 
namic configuration  of  the  axis,  even  In  forced  motion,  can  be  expressed  as 
a superposition 

4>ix,t)=^qx(t)\l>l{x)  [3] 

1 

/ where  <;.(  t )are  a set  of  generalized  coordinates.  For  these  motions  the  In- 
stantaneous kinetic  energy  of  the  body  can  be  given  as 


and  the  Instantaneous  potential  energy  as 


«,  Qi 


m 

lb) 
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!n  order  10  Illustrate  the  mode  functions,  we  may  c 1 u«  r t r »• 
special  case  of  a uniform  bar  of  length.  /.  Then  0#-lls  tu  m :«  function 
for  a rigid  vertical  translation,  04(jr)  - 1 - j dmotc;  tt.e  • < ■ 1 - . 

tatlon,  <*,(/)  denotes  the  mode  sha;  e when  the  bar  Is  flex'ng  *n  wo-i.oO.ed  v';r. 
lion,  t fie  center  moving  out  of  phase  with  the  two  enos,  etc.  .-.Is.  V„  uenote 
the  mass  of  the  bar,  Mx  denotes  the  moment  of  Inertia  of  the  bar  ..t . ut  ; 
mldlength,  etc. 


MOTION  OK  THE  WATEK 

Consider  now  the  motion  of  the  water.  We  assume  ti.at  tt.e  f.sw  1: 
incompress',  ve  and  tnv  Isold,  and  ti.at  a velocity  potential  #(  )<  x' s •:  s 


satisfies  the 

following  c 

ondlt Ions : 

-o 

at 

Interior  points  (,Q,C  of  tt.e  fluid 

.t'j 

*1 

1 

O 

at. 

Infinity 

17) 

e a " 

on 

tne  surface  of  the  bubble 

l$J 

where  a Is  an  outward  normal  and  Is  the  local  normal  velocity  of  tne  bubble 
outward;  and 

(«,1)  on  the  surface  of  tne  body 

I 

where  (n,A)1s  the  angle  between  the  vertical  and  the  outward  normal  to  tie  sur- 
face. Also  the  pressure  must  be  continuous  across  tht  bubble  surface  so  that, 
1 f Pj(Ols  the  Internal  pressure  In  tne  bubble  (assumed  uniform),  tl.en  at  its 
surface 

110) 

where  pk  's  the  hydrostatic  pressure  at  the  point  and  p Is  the  water  density. 
Equations  (o]  through  (10)  plus  an  equat'on  of  state  for  tne  gas  bubble  com- 
pletely specify  a boundary  value  problem  for  ♦.  This  function  may  now  be 
constructed  by  superposition  of  several  different  potentials  as  follows. 

First  tnere  Is  the  velocity  potential  ##  which  Is  due  to  the  motion 
of  the  gas  bubble  when  the  elastic  body  Is  replaced  by  an  Immovable  rigid 
surface  atul  the  bubble  Is  constrained  to  maintain  Its  original  motion.  Then 
Is  a harmonic  function  whose  gradient  vanishes  at  Infinity,  and 
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on  the  surface  of  the  bubble  G [11] 

on  the  surface  of  the  target  S [12] 

Second  there  are  the  velocity  potentials 

•i  ■ I -0.1,2  [13] 

Each  potential  describes  the  flow  when  the  body  Is  vibrating  In  a prescribed 
■ode  and  the  bubble  Is  absent.  Thus  #,  Is  a harmonic  function  whose  gradient 
vanishes  at  Infinity  and 

- coe  (a.  fc)  on  S [14] 

0 a 

With  these  definitions  it  is  easily  verified  that 

♦ MM 

I 

sat'sfles  Equations  lb]  through  [10]  except  perhaps  for  [8],  the  boundary 
condition  on  <»’.  In  order  to  satisfy  this  last  condition  it  would  be  necessary 
that 

- 0 on  (,’  lib] 

c n 

i condition  which  has  not  been  prescribed  In  order  to  avoid  making  0,  depend 

on  t.  h wevcr,  In  the  further  analysis  we  shall  suppose  that  the  bubble  Is 

aa 

sufficiently  far  from  the  body  that  the  values  of  at  (,’  are  small  and  — ' 

dn 

at  G Is  negligible. 

Now  we  can  write  the  kinetic  energy  of  the  water  as  an  Integral  over 
the  entire  volume  of  water  in  terms  of  the  potentials  and  #,  and  then  trans- 
form the  result,  by  means  of  Green's  theorem.  Into  integrals  over  the  sumfoee 
of  the  solid  .s’  and  the  surface  of  the  bubble  (». 


e n 


f*  - 0 

0 N 


T,  -SjWrfr- -/(#$:  a.  '£♦£> 

Substituting  for  ♦ from  [l1’]  and  using  the  boundary  conditions  ill], 
and  [14] 


M7] 

[12], 
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t 


0 cos  ( n.AI  d o * 

I 


P 

-» 


a a 


♦42' 


d<px 

I '<>  n 


)d<! 


Pi 


P 


S'  d*<\. 

- ‘i . a »*  Ma 


Some  of  these  terms  can  he  consolidated  by  means  of  a relation  deduced  from 
Green’s  second  theorem. 


)>< 

0v*>, 

b >1 

b '1 

- )> 

bi?, 

' fin 

*'  b it  1 

J, 

cos  (n.AMc 

1 $ ' 

tie  + 

I O c.ite 

1 1 

P *1 


i-  k 


which  may  he  substituted  In  the  first  term  In  11$]. 

Consider  now  t t.e  last  Integral  In  [ l $ ) . 
mately  uniform  overt.’  If  the  bubble  Is  sufficiently 
Also  averages  to  zero  over  t Hence  this  last 
compared,  for  example,  to  th.e  proceeding  term. 


The  factor  P,  is  approxl - 
distant  from  the  solid. 
Integral  can  be  neglected 


We  consider  now  the  Integrand  of  the  last  term  In  ].’0].  The  func- 
tion #,  describes  the  flow  which  would  prevail  If  the  bubo l e were  removed. 
Hence  Its  mean  value  ever  the  surface  of  (i  Is  equal  to  the  value  which  would 
exist  at  th.e  center  of  t However,  the  presence  of  the  burble  modifies  this 
flow  by  adding  a component  which  ’s  essentially  a dipole  term  proportional 
to  the  magnitude  of#,  . The  mean  value  of  this  Induced  component  over  t.’  will 
be  zero.  Also,  if  the  bubble  Is  sufficiently  distant  from  S.  th.e  bubble  must 
be  essentially  spherical,  and  r will  be  approximately  uniform  over  me  sur- 
face. Hence,  we  can  write  th.e  last  term  In  ICO]  as 


where  l’# 
exist  at 

| id)  can 


Is  the  volume  of  me  bubble  and  #(1,  is  th.e  value  for  w !o  would 
the  position  of  the  center  oft.’.  If  the  bubble  were  absent 
Likewise,  since  Is  almost  uniform  overt.’,  the  third  Integral  In. 
be  transformed 
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iL*,l’tda  ' T 


22] 


where  #w  Is  a suitable  mean  value  of  over  (.’.  Finally  we  can  further 
simplify  the  terms  In  [l8]  by  defining  an  entrained  mass  of  water 

9 ' r 

p|  0,  cos(M.U(ie  P ) 0,  co*(n, k)do 

The  symmetry  property  of  L,,  follows  from  Green's  first  equation. 

Substituting  (20],  [21  ] , [f.'],  and  [23]  Into  [18]  ind  rearranging 
terms,  the  kinetic  energy  of  the  water  now  becomes 


r,  ^ m,  + J +J't+  .1 


l 2*4  ] 


< .) 


Th.ere  ue  two  terms  which  contribute  to  the  potential  energy  of  the 
bubble  and  water.  First  there  Is  the  Internal  energy  change  of  tm  bubble 
In  Its  adiabatic  expansion.  This  depends  upon  the  volume  and  upon  the  equa- 
tion of  state  for  the  gas;  It  snail  be  denoted  simply  by  K(V  1.  second  there 

Is  the  potential  energy  of  the  water  whlc;  has  been  displaced  by  the  expanding 

bubble.  This  energy  must  equal  V’.  Thert  are  also  additional  potential 
energy  terms  which  are  due  to  the  migration  of  the  bubble  In  a field  In  which 

pressure  gradients  exist.  Such  gradients  will  be  due  to  gravity  and  to  thi 

presence  of  surfaces,  however,  these  additional  terms  will  be  neglected  * n 
order  to  simplify  the  analysis.  This  means  that  the  further  analysis  Is 
applicable  only  to  those  times  during  which  bubble  migration  negligible. 


Adding  all  terms  now,  we  have  for  tne  kltut'.v 
of  target,  water,  and  bubble, 


it  to. 


. t n 


1 il 


energte; 


T 


\ ' V, 


' 0 e 


? <M>  + o'Vv, 


• . ) 


* -J/  V,  ’ rt<  ^ 4 4 


. t 
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EQUATIONS  OF  MOTION  OF  THE  SUBMERGED  BODY 
We  can  now  ap^ly  Lagrange's  equation 

\ _U1  + ev  m 
dt  \dgj  a q,  dq, 

and  obtain  the  differential  equation  for  the  flexural  motion 

L tl  ij,  + P$,gVt+  M'-a'fq,-  = 0 
j 

(My+Lli)(q\  + (jfq,)  + Lxy  q]  = -p«h„  V„ 

where 


3*' 


M 


, , t _ » i 

' Mi  + Luai 


12  71 


12  8] 


[291 


It  appears  that  each  mode  of  motion  behaves  like  a linear  oscillator 

which  Is  Inertia-coupled  to  each  of  the  other  modes.  The  generalized  mass 

for  the  i**  mode  has  been  augmented  by  an  entrained  mass  of  water  La,  and 

this  Increased  mass  Is  also  effective  In  decreasing  the  apparent  natural 

frequency  according  to  Equation  [29].  The  generalized  force,  the  term  on  the 

• • 

right  In  Equation  [28],  Is  the  product  of  a time -dependent  factor  Vt,  which 
depends  only  on  the  bubble,  and  a space  factor  , which  depends  only  on 
the  flexural  properties  of  the  body  and  the  position  of  the  bubble.  An  alter- 
nate form  of  this  generalized  force  can  be  obtained  from  Equations  [20]  and 
[21],  whence  the  generalized  force  becomes 

Qi  = C08  (n,k)d«r  + I p<j,  da  [30] 

Jc  v dn 

si. 

Now  Is  approximately  uniform  over  G,  while  averages  to  zero  over  G. 

a . on 

Hence  the  last  term  Is  negligible.  Also  p+„  Is  the  pressure  due  to  the 

bubble  motion  when  the  target  Is  clamped  rigidly.  Hence  the  generalized 

force  can  be  computed  as  though  the  body  were  rigid. 

The  coefficients  L i}  ( 1 *y)measure  how  much  tne  Inertia  of  the  water 

couples  together  two  modes  which  had  been  Independent  for  vibrations  In  air. 

These  coefficients  may  vanish  If  there  Is  sufficient  symmetry.  For  example, 

suppose  the  target  Is  symmetric  about  a vertical  transverse  plane  through  the 

midlength.  Then  the  functions  0,  will  be  symmetrical  or  antlsymmetrical  about 

this  midplane  depending  on  whether  i Is  even  or  odd.  Likewise  0,  must  have 

the  same  symmetry  as  0< . Hence  for  such  a body,  whenever  i and  j are  of 

opposite  parity,  the  coupling  coefficient  L,j  must  vanish,  or 
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/•M  P o , tP , cos  ( u . k ) <1  o 0 


For  the  more  general  case  of  a body  of  arbitrary  shape  we  can  trans- 
form to  a new  set  of  coordinates  and  mode  shapes  so  that  In  these  new  coor- 
dinates there  w'll  be  no  coupling  coefficients.  Thus,  return  to  Equations 
and  lot  ] and  transform  to  a new  set  of  coordinates  (()  so  that  In 
terms  of  these  new  coordinates  the  kinetic  and  potential  energy  terms  which 
are  Independent  of  the  bubble  become  simultaneously  sums  of  squares. 


^ M,qi  + ^ L,,  q,  q,  - . V qr' 


'A  / , o» , v • ^ w ‘ o ‘ 

T 1 r 


o'nce  the  kinetic  energy  Is  a positive  definite  form,  this  can  be  accomplish- 
ed by  a real  l'near  transformation: 

-^a,rqr(t)  [51*] 


New  'f  we  define 


4>r  ' n.r^i 


*r  ^<J,r0,  [3o] 

I 

and  substitute  In  [2^1  and  [2o],  we  get  from  Lagrange's  equation 

1r  + u>/qr  P0rtVf 

where  Is  the  value  of  q> , at  the  position  of  the  center  of  G.  Th.e  details 
of  the  transformation  are  snown  In  the  Appendix.  Th.e  final  result  Is  that, 
each,  mode  of  motion  for  the  combined  system  of  water  and  body  acts  Indepen- 
dently as  a simple  linear  oscillator. 

Th.e  function  cannot  be  expressed  inalytlcally  as  a function  of 
the  bubble  poslt'on  because  of  the  generality  of  the  boundary  conditions. 
However,  we  can  derive  a useful  approximation  which  Is  applicable  over  a large 
range  of  particular  cases.  This  approximation  applies  to  submerged  bodies 
w:.en  (1)  th.e  distance  from  the  charge  to  a point  on  the  target  surface  Is 
large  compared  with  the  cross-sectional  dimensions  of  t he  target,  and  also. 
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(.  ) the  body  Is  symmetric  about  a horizon- 
tal plane  containing  Its  longitudinal  axis 
(s'nce  we  are  considering  vertical  motions). 

We  write  In  the  familiar  form 
of  equation  [38]  below,  In  terms  of  the 
values  of  and  d<p\/dn  at  S,  r„  denoting 
the  list  nee  from  (»  tc  . point  n S,  sec 
Figure  i . 


Figure  1 - Geometry  for 
Submerged  Body 


477-0 


IV 


-I, 


£ * 


[38] 


N-W  cut  the  body  Into  disks  of  width  dz  by  transverse  planes  perpendicular 
to  Its  axis  and  let  ds  be  an  element  of  the  perimeter  around  the  disk  at  r . 
Then,  since  — = 0,cos(HP). 


4 7T0ltf 


+ ) l i — cos  (n,k>  dsdx 

Jfl  J rg 


[39i 


Consider  the  special  case  for  which  rg  is  much  larger  than  the  crass  section 
of  the  disk,  and  let  r be  tr.e  distance  from  G’  to  the  axis  of  the  disk.  Then 


j a^T  (r^)  d*  =“i  0. 


l^o] 


or.  In  terms  of  the  unit  vectors  i , j , k , in  the  directions  of  the  (,*],( 
axis 


J*V  P f*,'  (n*j)ds  £ 141] 

Because  of  the  assumed  symmetry  about  the  horizontal  midplane,  ("n.i  ) and 

, •*>  ♦. 

\n.ji  must  both  be  symmetric  about  this  midplane,  while#,'  must  be  antisymme- 
tric about  the  same  section.  Hence  the  first  two  Integrals  must  vanish. 

Also  the  Integrand  of  the  third  term  In  [41 ] can  be  transformed  by  Equations 
[83]  and  [87]  of  the  Appendix,  so  that  [41]  reduces  to 

k £(£)<•-£  j (*;-*;) 

where  #,"ls  a particular  combination  of  the  functions  as  described  In  the 
Appendix.  Now,  returning  to  the  second  term  in  [39], 
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f -1  cos  Oi.lc)  dsi/jr  - I div  ^ dAdr=~-  ^ Adi  ^ ^ 

J 'tf  J ru  r' 

where  A la  the  section  area  of  the  disk  arid  the  left  side  may  be  considered 
is  an  Integral  over  the  entire  surface  of  the  disk.  The  second  Integral 
follows  from  Gauss'  theorem  applied  to  the  volume  of  the  disk.  Note  that  ( 

Is  the  depth  of  (1  below  the  axis  of  .S' , and  If  (i  Is  above  .S’,  then  £ Is  negative. 
Substituting  [ 4^2 J and  1 4 3 J into  [59J 


l n<p  ’ - 

iii 


ip‘  di 


AP  - in 


ip* ‘ dr 


144  j 


This  last  result  Is  a generalization  of  an  analysis  given  by  Lamb  for  the 
effect  at  a distance  due  to  the  rigid  translation  of  a solid  In  a liquid.*’ 
Finally,  substituting  this  value  for  </>  ,„  In  ( 37  J 


,1 

</,'+<*>]  1,’ 


4 Tt 


j"'  A />  - m 
l r* 


' d r 


145] 


Note  that  the  transformation  between  i/’and  q Is  the  Inverse  of  that  between  ip 
and  ip  as  detailed  In  the  Appendix.  In  many  practical  cases,  the  body  dis- 
places Its  own  weight  of  water,  so  that 

[(Ap-m)dx  *=  0 , 

Jo 

and  we  should  expect  the  second  Integral  on  the  right  to  be  small  compared 
with  the  first.  In  this  equation,  all  the  parameters  may  be  determined  ex- 
perimentally by  observing  the  free  vibrations  of  the  target  In  water  and  the 
gas  bubble  In  a free  field. 


1‘HOI’OKTIUNAL  BODY 

The  simplest  version  of  this  analysis  occurs  In  the  special  case 
where  the  density  of  the  body  Is  uniform  and  equal  to  p = w and  where  the 
body  Is  "proportional."  We  define  a proportional  body  as  one  for  which 

| <P,  cos  (n,  k)  dn  - a(  mip'  [4b] 


where  a,  Is  a constant,  Independent,  of  r.  Multiply  the  equation  by  ip,  and 
Integrate  over  the  length  of  the  body: 


I ip  <fi  <j>x  cos  tn.k)  dndr  = at  f mip  ip  dr 
■at  ’ •'  -'o  1 1 


Tir  rofor> 


//  : “ (l  P M 6 

t)  \r  I 1J 


147) 
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Hence  a proportional  Hotly  Is  one  for  which  the  cross  coupling  coefficients 
of  1 ) vanish,  and  the  transformation  to  new  coordinates  Is  trivial,  l.e., 


1,'  ~ <j,  <Mt  + !.  )* 

' * ‘ It 

/ .1/  \ i 

"•  (jr+V.  ) 


ii>  (M  + /,  r * \i> 

t i i • 


and  substituting  in  [ 4 5 J , 


«.-[iSr  1 ^ ,/r] i;' 


The  entrained  mass  /,  1(  has  disappeared  from  the  equation  except  for  its  effect 
on  the  natural  frequency.  This  notion  of  a proportional  body  Is  not  merely 
a mathematical  abstraction  but  It  Is  a useful  concept  In  practice  where  the 
mode  function  In  water  Is  experimentally  the  same  as  the  mode  function  In  air. 

A simple  consequence  of  this  applies  to  the  vertical  translation 
of  the  solid.  This  may  be  considered  a mode  for  which  i/>,  \J>a , a constant, 
and  w,,-*  0.  The  terms  in  [b8]  now  no  longer  have  any  dependence  upon  the 
elastic  properties  of  the  body,  and  they  may  equally  well  apply  to  a body  of 
water  with  an  envelope  Sot  the  same  shape  as  the  surface  of  the  body.  Hence 
the  vertical  motion  of  the  center  of  gravity  of  a proportional  body  is  the 
same  as  the  vertical  motion  of  the  centroid  of  the  displaced  water  In  the  ab- 
sence of  the  body. 

More  generally,  If  we  consider  a very  rigid  structure  of  arbitrary 
density,  then  all  modes  of  motion  should  be  negligible  except  those  which  de- 
scribe vertical  translation  and  rotation,  and  these  modes  cannot  be  appreciably 
coupled  to  the  flexure  modes.  Assume  sufficient  symmetry  to  eliminate  coupling 
between  translation  and  rotation,  and  denote  the  translation  mode  by  the  sub- 
script 0.  Then,  effectively  the  body  acts  as  a proportional  solid  for  this 
mode,  and  . , » | | 

* 1.  ‘ 


>1  “ </  hVI„  + /,„) 


Substituting  In  1 4 5 ] , we  get 


Considering  only  the  case  for  which  r»/ 


■K f 


w M,> 


J + /■„ 

4nrs  +1.„ 
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where  Mw  la  the  mass  of  the  displaced  water.  Since  everything  starts  from 
rest,  the  displacement  la  almply  proportional  to  the  bubble  volume.  This  sim- 
ple result  could  have  been  derived  much  more  directly,  without  any  reference 
to  elastic  motions,  by  baaing  It  on  the  equations  for  the  translation  of  a 
solid  In  a fluid.*' 


FLEXURAL  MOTIONS  OF  A FLOATING  BODY 

The  problem  of  the  flexural  reactions  of  a floating  body  to  a pul- 
sating bubble  can  be  treated  by  the  same  methods  described  for  the  submerged 
body.  An  added  complication  la  that,  besides  the  other  boundary  conditions, 
the  velocity  potential  ♦ must  vanish  at  the  tree  surface.  This  condition  Is 
necessary  In  order  to  approximately  satisfy  the  condition  that  the  pressure 
at  the  free  surface  must  always  equal  atmospheric  pressure. 

The  condition  that  — o at  all  points  on  the  free  surface  may 

be  satisfied  by  the  usual  device  of  constructing  an  Image  space  which  Is  the 
reflection  of  me  fluid  space  In  the  plane  of  the  free  surface,  and  where  the 
potential  at  any  point  Is  *he  negative  of  the  potential  at  the  corresponding 
point  In  the  fluid.  Hence  If  S Is  the  Image  of  the  underwater  surface  .S’  and 
(J  Is  the  Image  of  the  bubble  (!  (see  Figure  2),  then  and  must  satisfy 
the  additional  boundary  conditions 


Figure  - Geometry  for  a Flouting  Body 
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where 


i\,  on  0 - -Hg 


for  Image  point  on  G' 


and 


(a.  /cl  on  S - -oosUi.  A)  for  Image  point  on  S' 


l 54] 


[55) 


Then  the  same  analysis  as  before  leads  to  the  same  differential  equation, 
as  In  1371: 

<i,‘+  ul *</,'“  fi*,t  Vg-  x>J  ^ It>'  coabt.k)  do  l5o] 

and  i ‘ f . . , , . 

n,,  = +x>  i/>t  cos  (H.k)  do 

•\s 

Is  defined  as  an  Integral  over  the  underwater  area  only.  Likewise  the  gener- 
alised force  on  the  right  side  of  the  equation  Is  computed  by  Integrating  over 
>\  the  undeiwater  area,  only. 

However,  In  order  to  evaluate  0,^  by  Green's  function 


157  J 


the  integration  must  be  made  over  ,s'  and  ,s’.‘  This  leads  to  a factor  two  In  the 
equation  for  0^  under  conditions  analogous  to  [44J,  l.e., 


m 0 dx 


0L  (.*!/>-  m) 
r® 


d x 


15»J 


where  A Is  now  the  area  of  the  undei-water  section  of  the  body  at  a distance  x 
from  one  end. 

As  a spectal  example,  we  give  the  equation  for  the  vertical  flexural 
vibrations  of  a floating  proportional  body  of  uniform  density: 


«/,  + w,  1, 


2Ps  f‘ 

4 nM-  J# 


dx 

*1/ 


1 59  J 


Note  that  the  only  difference  as  compared  with  [48]  Is  the  factor  of  two  In 
t tie  generalized  force. 


MOTION  OF  THE  BUBBLE 

In  order  to  solve  the  differential  equations  for  the  motions  of  the 
body,  we  must  know  the  volume  acceleration  of  the  bubble  i,,  or  an  equivalent. 
In  a practical  case  this  might  be  known  from  prior  knowledge,  e.g.,  from 
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photographs  of  the  bubble  motion.  However,  the  equation  for  the  motion  of 
the  bubble  may  also  be  obtained  from  this  analysis. 

We  shall  assume  that  the  bubble  Is  spherical  at  all  times,  of  radius 
R(t),  and  that  near  the  bubble  the  velocity  potential  Is  given  by 


and 


160] 


161] 


Is  a harmonic  function  which  satisfies  the  boundary  condition  [11]  at  the 
bubble,  but  not  the  boundary  condtlon  [12 J at  the  body.  It  Is  clear  that  near 
the  body  the  potential  must  be  augmented  by  a term  representing  a "reflected" 
potential.  We  assume  that  the  body  Is  sufficiently  distant  from  the  bubble 
that  this  reflection  term  Is  negligible  at  the  bubble  Itself.  Then,  substi- 
tuting [6l]  Into  [2^]  and  1?6], 


' i.J 


ha 

2 


<?,  Q. 


+ 2ir  pR*R2  *Z  4 np^ 


R2R 


[62] 


V =£  T1"  “.2  «?  + | nIi*Ph  + E ( R ) 

I 

and  applying  Lagrange's  equations  to  the  coordinate  R, 


[65] 


£ 4rqiP*i#  R2  + 4npR3  R + 6rrp  R2  R2  + 4 irR2  ph  + |^-  - 0 [64] 


The  first  term  couples  the  pulsation  of  the  bubble  to  the  motion  of  the 
solid.  This  term  may  be  written  as  4 rrRtpS(1  where 

psc=£p*n>  mZp'<ii+i' 

X 1 

Is  the  pressure  at  G due  to  the  vibration  of  S.  Again  If  the  body  Is  suffi- 
ciently distant,  then  p5Cls  small  compared  for  example  with  pk  which  enters 
Into  the  fourth  term.  If  this  first  term  Is  neglected,  then  the  remainder 
of  Equation  [64]  Is  equivalent  to  the  differential  equation  which  Is  usually 
derived  for  a bubble  In  an  Infinite  fluid.'  The  equation  must  be  solved 
either  by  numerical  methods,  or  the  main  characteristics  of  the  pulsation 
can  be  estimated  by  methods  which  are  described  In  the  literature.3 
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EFFECT  OF  AN  INITIAL  SHOCK  WAVE 

The  foregoing  was  all  based  on  the  assumption  that  the  flow  which 
Is  generated  by  the  bubble  pulsation  Is  an  lncompresslve  motion.  This  Is  a 
common  assumption  In  bubble  theory  and  Is  known  to  be  accurate  during  a large 
part  of  the  bubble  oscillation  period  when  the  radial  velocity  of  the  bubble 
surface  Is  small  compared  with  t , the  velocity  of  sound  In  water.  However, 

In  the  Initial  stages  of  Its  motion,  when,  for  example,  a bubble  Is  created 
by  the  detonation  of  a charge,  a compressive  wave  Is  emitted  and  the  velocity 
potential  differs  markedly  from  a harmonic  function.  Vie  shall  assume  that 
this  flow  satisfies  the  wave  equation 


and  we  shall  determine  under  what  conditions  the  lncompresslve  analysis  re- 
mains valid.  It  Is  still  assumed  that  the  flows  described  by  ♦,  are  lncom- 
presslve. 

Lagrange's  equations  can  no  longer  be  applied  to  the  motion  of  the 
water,  because  It  has  been  demonstrated  only  that  they  are  applicable  when 
the  flow  Is  lncompresslve.3  We  return,  therefore,  to  the  mode  functions 
which  describe  the  normal  vibrations  In  air,  and  consider  the  forced  vibration 
In  the  i '*  mode  due  to  an  applied  pressure  of  magnitude 

f>4  - 0 + ^2,  » V,  l bo  ] 

J 

The  generalized  force  for  this  mode  must  be  computed  by  Integrating  this 
pressure  over  S with  - *l>,  coa  (w,  Hr)  as  a weighting  function.  Hence 

•W,  U,  + «/J  “ J 4>,  in.k)  do  + ^ p ^ j #,  ms  i H.k'do 

J l «>7] 

* p ft  £ % «■  ^ 1‘%)  ^ 

Again  we  use  Green's  second  Identity  to  transform  the  first  term  on  the  right 

’ I**  ^ do  “ J do  + | F2  <<t  [t>8] 


where  the  last  Integral  Is  taken  over  the  whole  volume  of  water.  Hence 


CONFIDENTIAL 


CONFIDENTIAL 


where  pU.H.C.t)  Is  the  local  excess  pressure  due  to  the  bubble.  Substltut’ng 

In  lt>7) 

•w.  + j PdT  17°) 

Now  Integrate  twice  with  respect  to  tine,  starting  with  time  zero  when  every- 
thing Is  at  rest  and  the  bubble  Is  created: 

M,  [l,  + a'  J Q.  dt  <*<]  + 7'..  <1.  P dr  [71] 

Consider  the  Integrand  of  the  last  term.  4,  Is  a function  which  drops  off 
as  from  the  target  (see  Equation  1^5)  )•  Also  It  Is  of  opposite  sign  above 
and  below  the  target.  Is  Initially  very  large  at  the  bubble,  but 

this  phase  travels  off  Into  the  water,  the  amplitude  decaying  as  -p  . 
Furthermore  It  Is  known  that  p becomes  negative  at  the  bubble  at  about  11 
percent  of  the  bubble  period  and  this  negative  phase  travels  out  Into  the 
water.  Hence  It  is  clea-  that  shortly  after  the  creation  of  the  bubble,  a 
t'me  I will  exist  at  which  the  last  term  In  Equation  [71]  will  be  negligible, 
or  even  zero,  compared  with  the  first  term  on  the  right,  and  after  this  time 
the  last  term  will  remain  negligible.  If  also  this  time  t is  small  compared 
with  the  natural  period  of  oscillation  In  the  i •*  mode,  i.e.,  if  a*f * « 4 jt*. 
then  the  Integral  term  on  the  left  side  of  this  equation  must  be  negligible 
compared  with  the  first  term.  Hence  at  this  time  t,  the  displacement  , as 
calculated  from  Equation  [71].  will  be  the  same  whether  or  not  the  compressl- 
billty  term  Is  included.  Furthermore,  after  this  time  the  compressibility 
effect  Is  known  to  be  negligible  and  need  not  be  Included. 

This  result  means  that  the  Incompress! ve  theory  Is  adequate  to  ex- 
p.aln  the  flexural  vibrations  of  the  body  even  though  the  solid  Is  Initially 
exposed  to  a compressive  shock  wave,  provided  that  the  period  of  oscillation 
is  long  compared  with  the  duration  of  the  compressive  wave.  On  the  other 
hand,  for  high-frequency  motions  of  the  body  the  generalized  forces  are  not 
related  to  Vt  In  a simple  way,  and  the  theory  Is  not  applicable. 

SUMMANY  AND  CONCLUSIONS 

We  may  sum  up  the  results  of  this  analysis  as  follows. 

1 . When  a flexible  body  is  submerged  In  water,  the  modes  of  flexure 
In  which  the  body  vibrates  become  modified  In  two  ways:  (a)  The  water 
increases  the  Inertia  of  the  system  and  so  decreases  the  natural  frequency, 
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and  (b)  the  Inertia  of  the  water  couples  together  the  various  modes  of  flexure. 
However  there  does  exist  a new  set  of  flexure  modes  In  each  of  which  the 
body  may  vibrate  Independently.  The  equations  for  these  new  modes  are  derived 
In  terms  of  the  old. 

2.  When  this  flexible  body  Is  exposed  to  a pulsating  bubble  In  the  li- 
quid, tne  bubble  motion  Is  affected  In  two  ways:  (a)  The  presence  of  the 
body  will  modify  the  pressure  field  from  the  bubble  and  thus  set  up  a pressure 
gradient  at  the  bubble  which  will  cause  the  bubble  to  migrate  to  the  body, 
and  (b)  the  mot  Ion  of  the  body  will  generate  a fluctuating  pressure  at  the 
bubble  which  will  modify  the  pulsation.  The  first  effect  has  been  Ignored; 
an  equation  has  been  derived  for  the  second  effect. 

).  The  bubble  will  generate  flexural  motions  In  the  body.  An  equation 
Is  derived  which  relates  the  motion  In  each  normal  mode  of  flexure  to  the 
volume  acceleration  of  the  bubble.  This  equation  Is  applicable  whenever  the 
bubble  Is  sufficiently  distant  from  the  body  so  that  (a)  the  normal  modes  of 
the  body  are  not  changed  by  the  presence  of  the  bubble,  (b)  the  migration  Is 
negligible,  and  (c)  the  distance  between  bubble  and  body  Is  large  compared 
with  the  section  diameter  of  the  body.  In  the  event  the  body  has  some  special 
symmetry  or  "proportional"  structure,  a simplified  form  of  these  equations 
Is  derived. 

4.  The  theory  Is  based  upon  the  assumption  that  the  flow  associated 
with  the  bubble  motion  Is  Incompresslve.  Yet  the  Initial  stage  of  the  motion 
may  Involve  a highly  compressive  wave  motion.  Nevertheless  It  Is  shown  that 
the  analysis  remains  applicable  provided  the  duration  of  the  compressive  wave 
Is  small  compared  with  (a)  the  period  of  pulsation  of  the  bubble  and  (b)  the 
period  of  vibration  of  the  body.  Another  effect  of  the  compressive  wave, 
which  has  not  been  discussed,  Is  that  It  provides  a mechanism  by  which  bubble 
energy  Is  radiated  away  from  the  system. 

It  Is  also  assumed  Implicitly  that  the  compressibility  of  the  water 
has  negligible  effect  upon  the  flexural  vibration  modes  of  the  body.  However, 
It  Is  clear  that  if  the  normal  velocities  of  the  surface  are  not  small  com- 
pared with  the  velocity  of  sound  In  water,  then  we  may  expect  that  the  effect 
of  the  water  on  the  normal  modes  Is  not  only  to  Increase  the  Inertia  of  the 
system  but  also  to  provide  a mechanism  for  energy  dissipation.  The  sound 
waves  radiated  by  the  vibration  will  drain  energy  from  the  bv.dy  and  thus 
dampen  the  vibration. 

The  experimental  verification  of  the  analysis  will  be  discussed  in 
a subsequent  report,  as  well  as  some  quantitative  implications  of  the  varia- 
tion of  flexure  with  target  and  charge  characteristics. 
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APPENDIX 

THE  RELATION*  BETWEEN  THE  EH EE  FLEXURAL  MODEL  IN  WATER  AND  IN  AIR 

We  derive  Here  aome  explicit  relations  between  the  normal  modes  and 
frequencies  of  a body  In  water  and  the  corresponding  modes  and  frequencies  In 
air.  These  relations  are  most  conveniently  obtained  In  matrix  notation.  We 
will  write  a matrix  In  bold  face  A,  Its  transpose  will  be  written  A uul  Its 
Inverse  as  A' . 

Let  q denote  the  column  matrix  with  elements  <j,(n. 

<>  denote  ttie  column  matrix  with  elements  tM.r), 

^denote  the  column  matrix  with  elements 

and  be  defined  only  on  surface  of  tr.e  body, 

M denote  the  diagonal  matrix  with  element.  M,, 

L denote  the  symmetric  matrix  with  elements  L,t  , ind 
a denote  the  diagonal  matrix  with  elements  a,. 

Then  from  Equations  [.’^J  and  [.'o]  we  can  write  for  the  kinetic  energy  and  the 
potential  energy  of  the  system  In  the  absence  of  the  bubble 

2T  - q (m  + l)  q = Q M a M '(m  + l)  o»  'M  <x  M d 1 72  j 


21'  q M o»  q (q  M a)  (a  M q) 


wnere  the  Insertions  have  been  made  In  order  to  express  V as  a sum  of  squares 
Now  diagonalize  the  matrix  In  brackets  by  means  of  an  orthogonal  1 ty  trans- 
format 1 on : 


M <*  '(m  + l)  a ’M  - ! R >A  R 


where  A Is  a matrix  such  that  A * Is  a diagonal  matrix  whose  elements  (w,  *) 

are  tne  characteristic  roots  of  the  matrix  In  brackets,  and  R Is  an  orthogon- 
al matrix  which  Is  compounded  of  the  characteristic  solutions  of  the  matrix 
In  brackets.  Hence 

2 T (d  M ' o*  R ')  A a (R  a M ■'  d)  - d' d'  1751 

2T  - (q  M e.R  'A ')A*  (A 'RaM-q)  q'A’q'  I7n] 

where  we  now  have  a new  set  of  normal  coordinates 


q'  [A  lR  o«  M ) q Aq 
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where  A Is  the  transformation  matrix  In  parenthesis.  Hence  V and  Tare  both 
now  sums  of  squares,  the  modes  are  Independent,  and  the  new  natural  frequen- 
cies are  the  Inverse  square  roots  of  the  characteristic  roots  of  the  matrix 
In  [74].  Also  we  define  an  associated  set  of  mode  functions  and  potentials 

l|>*  - A'l4>  • f-  A>  178] 


and  substituting  In  [72] 


where 


27*  * ^ (m  + l)  ^ q'M'<j'  + q*L'<|'  - V 


179) 


A 'MA 


*/ 

1 am  (b* 

.oT  T 


m it  x 


180] 


and 


and 


L*  « A 'LA  1 » +p)  (j)'  cos  (n.  A)  do 

X* 


[8l] 


M'  + L'  I the  unit  mutrlx. 


182] 


This  completes  the  specification  of  tne  new  coordinates  and  parameters  as 
unique  functions  of  the  old  coordinates  and  parameters. 

There  exists  a relation  between  Uie  new  potentials  and  the  new  mode 
functions  which  is  useful  In  the  analysis.  We  define  a new  series  of  functions 
0?lx)as  an  Integral  of  around  the  perimeter  » of  a section  of  the  body  at 
the  distance  x from  one  end,  l.e., 


( x)  cos  ( n , k)  ds 


[83] 


Let  denote  the  column  matrix  with  elements  and  consider  the  relation 
between  and  t|»'.  Now  since  the  functions  \l>,  , and  also  ’,  are  assumed  to 

can  be  expanded  in  terms  of  t/*,(x)  and 


form  a complete  set,  the  function 
we  can  take 


m ( 


--  m (,r)  B <]) ' 


184] 


where  B Is  an  undetermined  matrix  with  constant  coefficients.  Now  postmultiply 

both  of  above  by  the  matrix  <|>  * and  Integrate  with  respect  to  x. 

,*t_  - :i 

I <]>'  dx  B d>  4>*  oirfjr 
• 0 .’0 
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and  from  [8l  ] and  |d0]  and  substituting  In  l SU J 

F*  - 5 L*M'  'r  - (i-M')m'  *♦' 

$'  - p (♦**  - r) 


where 


>"  - M*  l^' 


and  represents  a column  matrix  of  mode  functions  which  are  particular  combina- 
tions of  the  mode  functions  This  final  result  Is  an  expression  for  the 

potential  distribution  on  the  surface  of  the  body  as  a function  of  the  motions 
which  generate  the  potential.  It  Is  Interesting  that,  In  general,  the  poten- 
tial ^(‘  does  not  depend  on  the  single  mode  but  also  on  all  the  modes  \1>S\ 
j M . This  dependence  on  all  the  modes  enters  through  the  presence  of  tit," 

In  18?]. 
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